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Abstract 

In this paper we study a class of stochastic partial differential equations in the whole space 
R d , with arbitrary dimension d > 1, driven by a Gaussian noise white in time and correlated 
in space. The differential operator is a fractional derivative operator. We show the existence, 
uniqueness and Holder's regularity of the solution. Then by means of Malliavin calculus, we 
prove that the law of the solution has a smooth density with respect to the Lebesgue measure. 
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1 Introduction and general framework 



In this paper we consider the following stochastic partial differential equation (SPDE for abbrevi- 
ation) given by 

dii 

— (t, x) = Vfu(t, x) + b(u(t, x)) + a(u(t, x))F{t, x), ^ 
u(0,x) = uo(x), 

where (t,x) € [0, T] x R d , d > 1, a = (a%, . . . , a^), 5 = (Si, ...,5^) and Vf denotes a non-local 
fractional differential operator to be defined below. More precisely, we will assume along this paper 
that 

on e]0, 2] \ {1} and |<5j| < min{aj, 2 — cti}, i = 1, . . . , d. 

The noise F(t,x) is a martingale measure (in the sense given by Walsh in |3U| ) to be defined with 
more details in the sequel. The coefficients b and a : R — > R are given functions. We shall refer to 
Eq. (Hm) as Eqf(d,b,(r). 

Here, we are interested in solutions which are real valued stochastic processes. Firstly, we 
establish sufficient conditions ensuring the existence, the uniqueness and the Holder regularity of 
solutions u(t, x) whenever they exist. Secondly, we study the existence and smoothness of the 
density of the law of the solution u(t, x) for fixed t > and x G R rf . 

Our results of existence, uniqueness and regularity of the solution extend those obtained in [2j for 
Eqf(l,b,a) and P for Eq^ (1, b, a) both driven by a space-time white noise W and 1 < a < 2. 
However in their framework, the regularity of the noise rises up a > 1 as a necessary condition 
for the existence of the stochastic integral as an L 2 (Q) random variable. This problem does not 
appear in our case (a G]0, 2] \ {1}) because the noise F is smoother than W . 

The equation Eq2(d,b,a) recovers for instance the stochastic heat equation in spatial dimension 
d > 1 studied by many authors among others jH], [IS]. With the notations adopted above it corre- 
sponds to Eq$(d, 6, a). Indeed, when 5, = and ai = 2 for i = 1, . . . , d the operator Vq coincides 
with the classical Laplacian operator in R rf . 

Various physical phenomena involving diffusion and interaction of particles can be described by 
the equations Eqf(d,b, a) when suitable assumptions are made on the coefficient b and a. The 
non-local property in this equation is due to the presence of T>f and the non-linearity comes from 
the general form of b and a. This equations can also be interpreted as random perturbation of 
deterministic equations Eqf(d,b,0) by multiplicative noise of the form a(u(t, x))F (t, x) . 
In probabilistic terms, replacing the Laplacian by its fractional power (which is an integro-differential 
operator) leads to interesting and largely open questions of extensions of results for Brownian mo- 
tion driven stochastic equations to those driven by Levy stable processes. In the physical literature, 
such fractal anomalous diffusions have been recently enthusiastically embraced by a slew of inves- 
tigators in the context of hydrodynamics, acoustics, trapping effects in surface diffusion, statistical 
mechanics, relaxation phenomena, and biology (see e.g. |2H|) EHj- [HI]] [111; ffil)- 
A probabilistic approach to the equation EqQ(d,b,0) is made by means of the Feynman-Kac for- 
mula (see 0, [1], J2])- Solutions to other particular fractional differential equations are given as 
functionals of stable subordinators. This representation provides explicit form to the density of the 
3/2-stable law and to the density of escaping island vicinity in vortex medium, see [S]. In other 
words the laws of stable Levy processes or stable subordinators satisfies fractional equations like 
E<$(d,b,0). 

SPDE lies at the intersection of several different areas, and different groups have taken the 
subject in their own directions. For example, one could start with the functional analytic approach 
to PDE, and develop a similar framework for SPDE. Reasoning this way, one views solutions of 
SPDE as random variables taking values in Sobolev spaces (or other spaces), and the terms in the 
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equations as operators and operator-valued random variables. 

Another possibility is to think of an SPDE as a stochastic process. In this approach, we focus on 
the intriguing properties of Gaussian processes and martingale measures. We will call this point of 
view the Walsh approach. 

To give a flavor of this approach, let us consider two particular processes, the first one is the 
superprocess or Dawson-Watanabe process. Roughly speaking, the density of the its law of this 
solves the SPDE 

-^(t, x) = T%u(t, x) + y/u(t,x)F(t, x) (d = 1, Vl = A) 
and the superprocess considered in ^1 whose associated SPDE is given by 

^(t,x) = Lu(t,x) - ^^-(t,x)W t + y/u(t,x)F(t,x) 

where L is a second order differential operator, g is a given function and W is a Brownian motion. 
Then super-processes arise as a limiting process in population dynamics and genetics. This con- 
nection leads to an analysis of certain SPDE in terms of particles. One of the sources of SPDE was 
the Zakai equation of filtering theory, and there was some work on the stochastic Navier-Stokes 
equation, but the theory of SPDE itself was not closely tied to applications. Luckily, physicists and 
engineers have now caught up, and SPDE is appearing in more and more scientific fields. We may 
mention some of these connections, such as randomly growing surfaces in solid-state physics and 
the theory of randomly moving polymers. Also, the stochastic Navier-Stokes equation has received 
much more intensive study in recent years. 

Our own motivation was, however, the connection between fractional operators (fractional power 
of the Laplacian) and super-processes. For example equation of the form Eq$(d, b, a) are related to 
some measure-valued processes. More precisely Eq^l, 0, is satisfied by the density for the law 
of a measure- valued branching diffusion for 1 < a < 2 and 5 = 0. Another example comes from 
Fleming- Viot models, that is, under some conditions on this diffusion process its law is absolutely 
continuous w.r.t. the Lebesgue measure with continuous density u(t, x) for each (t, x) G]0, +oo[xR. 
Moreover it satisfies the following SPDE 

— it, x) = V%u(t, x) + y/u(t, x)F{t, x) - u(t, x) / y/u(t,y)F(t,y)dy, 

at JR 

for d = 1, 1 < a < 2 and F(t, x) is a space-time white noise in x M. For complete description of 
measure-valued branching and Fleming- Viot diffusions the reader may consult |14j (see also and 
The above equation is the perturbation of Eqfi^l, 0, \f~-) by u(t, x) f R (u(t, y)) 1 / 2 F(t, y)dy. To 
recover this equation in Eqf(d, b, a) we allow a to depend on the whole trajectory of u. Nevertheless 
we are not studying this problem in this paper and restrict our selves to Lipschitz coefficients. 
The case where < a < 1 was motivated by the work [3] where he attempts to use the 2/3- 
order fractional Laplacian modeling of enhanced diffusing movements of random turbulent particle 
resulting from non-linear inertial interactions. A combined effect of the inertial interactions and the 
molecule Brownian diffusivities is found to be the bi-fractal mechanism behind multifractal scaling 
in the inertial range of scales of moderate Reynolds number turbulence. Accordingly, a stochastic 
equation is proposed to describe turbulence intermittency. The 2/3-order fractional Laplacian 
representation is also used to construct a fractional Reynolds equation for non-linear interactions 
of fluctuating velocity components, underlying turbulence space-time fractal structures of Levy 
2/3-stable distribution. The new perspective of his study is that the fractional calculus is an 
effective approach modeling of chaotic fractal phenomena induced by non-linear interactions. 
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Fractional derivatives and integrals, usually known as fractional calculus, have many uses and 
they themselves have arisen from certain requirements in applications, such as fractional integro- 
differentiation which has now become a significant topic in mathematical analysis. It has applica- 
tions in various fields namely quantitative biology, electrochemistry, transport theory, probability 
and potential theory to mention a few. We refer the reader for a complete survey on the fractional 
integrals and derivatives to |24| and [2U] (and the references therein). In recent years, fractional 
derivative has long been found to be a very effective means to describe the anomalous attenuation 
behaviors, some parametric seismic wave propagation have been model via the fractional deriva- 
tive model. Here will not discus this topic and limit our selves to the above connections between 
fractional calculus and mathematical and physical phenomena. 

In the literature, various fractional differential operators are defined (see ^U], ^H]) ED J- The re- 
sults of this paper apply to several of them, such as fractional Laplacian, Nishimoto operator and 
the non-self adjoint fractional operator introduced in 13. and used in |2] to the study of stochastic 
partial differential equation. 

Basic notations, definitions and preliminary results of the operator T>f and the noise F will be 
presented in the following subsequent two subsections. The section 2 is devoted to the existence 
and uniqueness and the section 3 deals with the Holder regularity result. The section 4 contains 
the existence and smoothness of the density of the law of the solution. In appendix we prove some 
technical results which will be used in the proofs. 

The value of the constants along this article may change from line to line and some of the 
standing parameters are not always indicated. 

1.1 The operator U$ 

For the sake of notational simplicity we denote by D? the fractional differential operator in dimen- 
sion d = 1 and T>f in dimension al > 2. 

Let us give the definition of the operator D? (i.e. in one space dimension). 

Definition 1.1 The fractional differential operator Df is a non-local operator defined via its 
Fourier transform T by 

F(D?<p) (0 = — lei" exp (-«5|sgn(0) ^(v) (0, 

where i 2 + 1 = 0. 

The operator Dg is a closed, densely defined operator on L 2 (M) and it is the infinitesimal 
generator of a semigroup which is in general not symmetric and not a contraction. It is self adjoint 
only when 5 = and in this case, it coincides with the fractional power of the Laplacian. 

According to 13 , Df can be represented for 1 < a < 2, by 

DMx) _ r- + -yw-ww Uti { .^ } (y) + 4 W.M) * 
J —oo \y\ 

and for < a < 1, by 

DMx) = f°° * {x ± y .l~ Q ^ {x) U s _i hoo , 0) (y) + 4i ( o, + oo)(y)) dy 
j —oo \y\ 

where «£. and k+ are two non-negative constants satisfying k s __ + > and tp is a smooth 
function for which the integral exists, and (p' is its derivative. This representation identifies it as 
the infinitesimal generator for a non-symmetric a-stable Levy process. 
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Let G a> s(t,x) denotes the fundamental solution of the equation Eqf (1,0,0) that is, the unique 
solution of the Cauchy problem 

( du 

J Qt(t,x) = Dfu(t,x), 

\ u(0,x) = 8 Q (x), t > 0, x G R, 

where 6q is the Dirac distribution. Using Fourier's calculus we obtain 

G a ,s(t,x) = J exp (-izx - t|z| a exp (-z<5^-sgn(z))) dz. 

The relevant parameters, a, called the index of stability and S (related to the asymmetry) improp- 
erly referred to as the skewness are real numbers satisfying a g]0, 2] and \S\ < {a, 2 — a}. 
The function G aj( $(i, •) has the following properties. 

Lemma 1.1 For a G (0, 2] \ {1} such that \S\ < min(a, 2 — a) 

(i) The function G Uj s{t, •) is the density of a Levy a-stable process in time t. 

(ii) The function G aj g(t,x) is not in general symmetric relatively to x. 

(Hi) Semigroup property: G aj $(t,x) satisfies the Chapman Kolmogorov equation, i.e. for < s < t 

/+oo 
G a ,s(^y) G a,s( s ' x - y) d v- 
-oo 

(iv) Scaling property: G a> s(t, x) = i -1 / a G Qj< $(l, t~ l l a x). 

(v) There exists a constant c a such that < G Qj( $(l,a:) < c Q (l + |x| 1+a ), for all 

For d > 1 and any multi index a = (a±, . . . , ay) and 5 = (S±, . . . ,Sd), define the operator V$ by 

i=i 

where D?' denotes the fractional differential derivative w.r.t the i—th coordinate. 

Let G a> s(t,x) be the Green function of the deterministic equation Eqf(d, 0,0). Clearly 

d 

G ay s(t,x) = JjG aii<5i (t,Xi) 
i=i 

= J^jd f Rd e *V ^-*(?' x ) -*^l&Texp (-^|sgn(&))^ df, 
where (•, ■) stands for the inner product in R d . 
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1.2 The driving noise F 



Let <S(M + ) be the space of Schwartz test functions. On a complete probability space Q,P), 
the noise F = {F(ip), ip G is assumed to be an L 2 (Q, G, P)-valued Gaussian process with 

mean zero and covariance functional given by 

J(<p,i/})= ds T(dx) (<p(s,-) * $(s,-)) (x), (p,tpeS 

Jr + JR d v 7 

where ip(s,x) = ip(s,—x) and T is a non-negative and non-negative definite tempered measure, 
therefore symmetric. Let [i denote the spectral measure of T, which is also a tempered measure 
(see that is \i = J r_1 (r) and this gives 



J{<p,fJ)) = ds /z(de)^(s,-)(0^(*»-)(0, (1-2) 

JR+ JR d 

where z is the complex conjugate of z. 

Following the same approach in [S] , the Gaussian process F can be extended to a worthy martingale 
measure M = {M(t, A) := F([0,t] x A) : t > 0, A G B b (M. d )} which shall acts as integrator, in 
the sense of Walsh [HOI, where i3b(M d ) denotes the bounded Borel subsets of M d . Let Qt be the 
completion of the u-field generated by the random variables {M(s,A), < s < t, A G Bb(M. d )}. 
The properties of F ensure that the process M = {M(t, A), t > 0, A G Bb(R d )}, is a martingale 
with respect to the filtration {Q t : t > 0}. 

Then one can give a rigorous meaning to solution of equation Eqf(d, b, a), by means of a jointly 
measurable and ^-adapted process {u(t, x) : (t, x) G M + x M. d } satisfying, for each t > 0, a.s. for 
almost x G the following evolution equation: 



u(t,x) = / G at s(t,x - y)u (y)dy 

JR d 

+ ds G a ,s(t- s,x - y)b(u(s,y))dy 
JO JR d 

+ 11 G atS (t-s,x-y)a(u{s,y))M{ds,dy) (1.3) 
JO JR d 

Throughout this paper we adopt the following definition. 

Definition 1.2 A stochastic process u defined on Q x M.+ x M d , which is jointly measurable and 
Qt~adapted, is said to be a solution to the fractional SPDE Eq?(d, b, a), if it is an M-valued fields 
which satisfies (j l.ty and 

sup sup E\u(t, x)\ p < +oo, for some p > 2. 

t€[0,T] x€R d 

Let us formulate our assumption concerning the fractional differential operator Df and the 
correlation of the noise M. 

If we take a = 1 and use the formula (|1.2fl the stochastic integral 

f [ G at s(T-s,x-y)M(ds,dy) 
Jo JR d 

appearing in (|1.3|) is well defined if and only if 



f T ds [ v(dZ)\FG a , s (s,-)(0\ 2 < +oc. 

Jo JR d 
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Indeed 



/ / G a)S (t-s,x-y)M(ds,dy) 
Jo JR d 

ds [ r(<fe)(G a ,i(t - s, •) * G^(t - a, - *) 
f ds f ti(dti)\FG a>5 (t-s,-)ti)\\ (^ = T). 



For a given multi index a = (a\, . . . , aj) such that < aj < 2, i = 1, . . . , d and any £ G M d we use 
the notation 



s«(o=i>r- 



i=i 



The following lemma gives sufficient condition for the existence of stochastic integrals w.r.t. the 
martingale measure M. 

Lemma 1.2 There exist two positive constants c\ and c 2 such that 



ci 



< 



1 + s a (0 

Proof. By the definition of G a g we have 



f [ nm\FG a ,s(s,-)(0\ 2 ds<c 2 [ ^f] . (1.4) 



Therefore, 



FG a , s (t, •)(£) = exp ^-i^l^r exp (-t^agu^))! 
jf*|^G a) i(a,-)(0| 2 cto = *f2t53|Ci| ai cos , 



(1.5) 



where $(x) = i(l — e :c )/x. Since $ is a decreasing function, we get 



*(2tS a (£)) < $ ^2t^|^|^cos(^|)^ <*(2tAs5 a (0) 



(1.6) 



with k = mini<j< ( i(cos((5j7r/2)) > since — ir/2 < 5iir/2 < ir/2. The mean value theorem applied 
to the function (1 + x)(l — e~ x ) yields 

t t / \ 2t 
< *(s) < 



1 + x y ' 1 + x 

for all x > and f G [0, T]. This together with (fT3|) and lfO|l imply that for all t € [0, T] 



1 + 2*S Q (£) 



< f \FG a , s (s,-)(0\ 2 ds< - 
Jo 1 



2f 



+ 2tK5 a (0" 



The conclusion follows by taking the integral over M. d with respect to the measure \x. 



□ 
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Remark 1.1 The upper and lower bounds in fli.^D do not depend on the parameter 5. When a, = 2 
for all i = 1, . . . , d then ^(C) = ||£|| 2 an d the bounds in \l-4\ are the same ones which appeared in 
0/ (see also \25}/). that is 



i + IKII 



< +00. 



Our main assumption on the noise relies an the integrability condition w.r.t. the spectral 
measure (i. 
Assumption (H") 



/ n ^fl^v < +°° where 77 €]0,1]. 

J R d (1 + £> a [€)) 



This condition states that there are not too many high frequencies in the noise, which turns out to 
be reformulated into a condition on the integrability of T (measure characterizing the covariance 
of the noise), near the origin in the case where T is a non-negative measure on M. d . The condition 
(Hg ) means that [i is a finite measure, which is equivalent to say that T is a uniformly continuous 
and bounded function on M. d . For this reason we do not consider this case in this paper. 
When /i is the Lebesgue measure on W d which is the spectral measure of the white noise on W d that is 

the noise corresponding to T = 5q, the condition (H^) is equivalent to f Q 00 r 1=1 Q * (1 + r) v dr < 

+00 which is finite in and only if r\ > X^=i(V a «)- 

1.2.1 Some examples 

Let us now, consider a class of covariances measures T for which the condition (H^) gives an optimal 
criterion. In this part we assume cti = 2, % = 1, . . . , d and use \x\ as the Euclidian norm of x. 



1. Riesz kernels. Let T(dx) = \x\ ^dx, 7 e]0, <i[. It is known that the spectral measure is fi(dx) 



c|x| 7 d dx. Then (H 2 ) is satisfied if and only if 7 e]0, 2r] A d[. 



2. Assume that the spectral measure /U is either finite or absolutely continuous with respect to 

Lebesgue measure and \i{dx)jdx = p(x) £ L p (M. d ) for some p > 1. For \i finite (H 2 ) holds 
for any 7 g]0, 1[. In the second case, (H 2 ) holds if d{p -l)/2p < tj < 1. That is if p £ 
[l,d/(d-2)Vl]. 

3. Brownian free field. Let J(ip,ip) = ((—A + m 2 )" 1 ^, V>) , where A is the Laplace operator in M. d 

and m > 0. In this case fi(dx) = (2ir)~ d / 2 (\x\ 2 + m 2 )~ l dx. Therefore, (H 2 ) holds if and only 
if 77 £](d - 2)+/2 A 1, 1[. That is, for d = 1, 2, (H 2 ) holds for any 77 € (0, 1), for d = 3, (H 2 ) 
holds for any n G]l/2, 1[; for d > 4, (H 2 ) does not hold. 

4. Bessel Kernel. Assume > and let T be the Bessel kernel of order 0. That is T(dx) = Bp{x)dx 

where Bp{x) = vp^ J +o ° t~ 1Jr ^~^f 2 e~ t ~^l u dt, x eM. d and is a constant. The spectral 
measure /j, of T is given by r(Bp)(£) = (1 + |£| 2 ) -/3/2 . Therefore (H 2 ) is satisfied if and only 
if 

C c+00 

/ (1 + lei 2 )-"-^ 2 ^ = c d (1 + t)-^i 2 t- l+d ' 2 dt 

JVL d Jo 

which is converg ent for d<2r] + f3. Hence (H 2 ) holds for n £}(d - /3) + /2, 1[ 
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2 Existence and uniqueness of the solution 

We shall also need a hypothesis on the initial condition uq: 

(H.l) no is an C?o _m easurable random variable s.t. sup^^d E\uq(x)\ Po < +00 for some po large 
enough in [2, +00 [. 

The main result of this section is the following: 

Theorem 2.1 Suppose that b and a are Lipschitz functions and assume that (H.l) and (H") are 
satisfied. Then, there exists a unique solution u(t, x) to Qi.ffi such that 

sup sup E\u(t,x)\ p < +00, for any T > and p E [2,Po\- (2-7) 

t£{0,T}x€R d 

Remark 2.1 Our results improve those of J^j and ^jj to d-dimensional case and extend the result 
in 1$ to the fractional setting, both with on in ]0, 2] \ {1} for i = 1, . . . , d. 

Proof. The proof of the existence is done by the Picard iteration scheme. That is, we define 
recursively for all n > 1. 

u°(t,x)=[ G a ,s(t,x - y)u {y)dy, 

t 

< u n (t,x) =u°(t,x)+ / G ajS (t- s,x-y)a(u n ~ l {s,y))M(ds,dy) (2.8) 
t Jo JR d 

+ ds G Qi s(t - s,x - y)b (u n-1 (s,2/)) dy. 
\ JO JR d 

Fix p £ [2,po\. To prove the theorem we only need the following facts: 

(i) For n > 0, u n = {u n (t,x) : (t,x) G [0,7] x R d } is predictable process which satisfies: 

(P n ) sup sup E\u n (t,x)\ 2 < +00. 
te[o,T] x m d 

This proves that the sequence u n is well defined. 

(ii) For T > and any p G [2,p ], sup„> sup t6[0)T] sup^d E\u n (t,x)\ p < +00. 

(iii) Set v n (t) = sup xm d E \u n+1 {t,x) -u n (t,x)\ p , n > 0. Then 

v n (t)<cf v n _ 1 {s){J{t-s) + l)ds, (2.9) 
Jo 

where 

J(t)= f \FG a>5 (t,-)(0\ 2 KdO- 

JR d 

From this we conclude that the sequence u n (t,x) converges in L p (£l). 
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Statement (i) is checked by induction on n as a consequence of the hypotheses (H.l), (H") and the 
properties on the coefficients. Indeed Lemma lA.ll applied to the probability measure G a> $(t, x—y)dy 
gives 

2 



£|u u (i,x) 



E 



G a> s(t,x - y)u (y)dy 

< E G a> s(t,x - y)\u (y)\ 2 dy 

< sup E\uo{y)\ 2 < sup E\uo(y)\ Pa < +00. 



This proves (Pq). Now, we assume that the property (jPfc) is true for any integer t < n — 1 (n > 2). 
We want to check (P n ). Burkholder's inequality, the induction hypothesis the Lipschitz condition 
on a and (H") imply 



sup sup E 

t£[0,T] xi 



G a>s (t -s,x- y)a [u n ~\s, y)) M{ds, dy) 



< c sup snpE(l + \u n 1 (s,y)\ 

sS[0,T] y<ER d 



l ) f da [ \TG a , s (s,-)(0\ 2 KdO- 
' Jo Jm d 



For the term containing the drift coefficient b we apply Lemma lA.ll with respect to the measure 
G at s(s,y)dsdy to obtain 



sup sup E 

te[o,T] x m d 



ds G ajS (t - s, x - y)b (u n (s,y)) dy 



-;[0,T] 

which is also finite. Hence 



< c sup sup E (l + \u n 1 (t,x)| 2 ) ( / ds G a) s(s,y)dy 
te\o,T]xm d ^ ' \h JR d 



sup sup E\u n (t,x)\ 2 < c sup sup£Ml + |u n 1 (t,x)\"), 
te[o,T] x m d te[o,T] x m d ^ 



we deduce (P n )- 

It is easy to check that u n is predictable. 
Statement (ii): Fix p £ [2,po\. We first prove that for any n > 1, t € [0, T] 

sup E\u n (t,x)\ p < a +c 2 f sup £ |ii n " 1 (s,x)| P (J(t - s) + l)da. 

The arguments are not very far from those used in the proof of (i). Indeed, we have 

E\u n (t,x)\ p < c(C (t,x) + A n (t,x) + B n (t,x)) 

with 

C (t,x) = E\u°(t,x)\ p , 



(2.10) 



(2.11) 



A n (t, x) 
B n (t,x) 



E 
E 



G atS {t -s,x- y)a (u n \s, y)) M(ds, dy) 



ds G ajS (t - s, x - y)b (u 11 1 (s,y))dy 

JO JR d 
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(2.13) 



Jensen's inequality and assumption (H.l) yield 

sup sup C (t,x)< sup sup (E\u°(t,x)\ po ) p ° < +00. (2.12) 

t£[o,T] x m d te[o,T]xm d 

Burkholder's inequality and the linear growth condition of a lead to 

sup A n (t,x) < c I ( 1 + sup E \u n ~ l {s,y)\ P J J(t - s)ds. 
Moreover, Lemma |A. II the linear growth condition of b and (i) of Lemma ll. II imply 
sup B n (t,x) < c ds 1 + sup E |ii n_1 (s,y)| P / G ajS (t - s,y)dy 

x£R d JO \ y£R d / J~R d 

= cf 1 1+ sup E\u n - 1 (s,y)\ P ) ds. 
Jo \ y m d J 

Plugging the estimates (|2~T2l to (ETUI) into (EHTj) yields (0TTT71) . 

Finally, the conclusion of part (ii) follows by applying Lemma |A . 21 a uoted in the appendix below to 
the situation: f n (t) = sup x< z R d E\u n (t,x)\ p , k\ = ci, k,2 = 0, g(s) = C2(J(s) + 1), with c\,C2 given 
in (I2~TI7|) . 

Statement (hi): Consider the decomposition 



(2.14) 



E\u n+1 (t,x) -u n (t,x)\ p < c(a n (t,x) + b n (t,x)), 



with 



a n (t,x) 



E 



Jm. d 
t 

ds 







G a , s (t -s,x-y) [a(u n (s, y)) - a{u n -\s, y))} M (ds, dy) 
G a , s (t -s,x-y) [b(u n (s, y)) - b(u n -\s, y))] 



b n (t,x) = E 

Again Burkholder's inequality, Lemma I A . 1 1 and the Lipschitz property of b and a imply that 



E\u n+1 (t,x) -u n (t,x)r < c 



f ( sup E\u n {s,y)-u n - l {s,y)\ P \ ( J(t - s) + 1) ds. 
Jo \ y eR d J 



This yields (|2.9|) . We finish the proof by applying Lemma rA.2l in the particular case: f n {t) = v n (t), 
k± = &2 = 0, g(s) = c(J(s) + 1), with c given in (|2.9j) . Notice that the results proved in part 
(ii) show that v := sup 0<s<T fo(s) is finite. Hence the series Sn>o( u ™(^)) 1 ^ P converges for any 
p G [2,po]- We then conclude that {u n (t, x) : (t,x) G [0,T] x M. d } converges uniformly in L P (Q) to 
a limit u = {u(t, x) : (t,x) G [0,T] x R d }. It is not difficult to check that u satisfies conditions of 
the Definition 11.21 and therefore the existence is completely proved. □ 

Let us now give a short proof for the uniqueness. 

Let u\ and U2 be two solutions to (|1.3|) 

F(t,x) := E\ui{t,x) - u 2 (t,x)\ 2 

and H(t) = sup zgIK d F(t, x). Then the isometry property for stochastic integrals, Lemma [A. II and 
the Lipschitz condition on b and a yield 

H[t) < c [ H(s)(J(t -s) + l)ds. 
Jo 

By iterating this formula and using Lemma [A. 21 we deduce that H = 0, hence, ui(t,x) = U2(t,x) 
t, x a.e. □ 
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3 Path regularity of the solution 

In this section we analyze the path Holder regularity of u(t,x). The next Theorem extends and 
improves similar results known for stochastic heat equation, corresponding to the case an = 2, 
<5j = 0, i = 1, . . . , d. (see for instance |26j^. Contrary to the factorization method usually used in 
high dimension (see e.g. and [7j) we use a direct method to prove our regularity results; in 
which the Fourier transform and the representation of the Green kernel play a crucial role. 

Let 7 = (71,72) such that 71,72 > and let K be a compact set of M. d . We denote by 
C 7 ([0, T] x K;M.) the set of 7-Holder continuous functions equipped with the norm defined by: 

11/11 \tu m 1 \f(t,x) - f(s,y)\ 

\\f\\ 7 ,T,K = SU P l/MI + SU P SU P tt -m ~ „ — -jiw 

(t,x)e[o,T]xK s^te[o,T]xjtyeK l J - s \ +\\x-y\\ 

(H.2) There exists p G]0, 1[ s.t. for all z € K (compact subset of W 1 ) 

sup E\u (y + z)- u (y)\ P0 < c\z\ PPo 



for some po large enough in [2, where ao = mini^j^^laj}. 

The main result of this section is 

Theorem 3.1 Suppose that b and a are Lipschitz functions. Assume that (H.l), (H.2) and (H^) 
hold for 7] e]0, 1[ and let u be a solution to equation Then u belongs a.s. to the Holder space 

C 7 ([0,T] xir;M) /or 0< 71 <min{^ti(pM),(l-??)/2}, 0< 72 <min{ /0 ,a (l-??)/2, 1/2} and 
K compact subset ofM. d . Moreover E ||"u||^ T ^- < +00 for any p £ [2,po]- 

In what follows, we establish some technical and useful results that will be needed in the proof 
of the regularity of the solution. 

Let uq be a given random fields satisfying (H.l). Set 

u°(t,x)= G a j(t,x - y)u (y)dy. 
Proposition 3.1 Suppose that uo satisfies (H.2) then 

(i) For each x G R d a.s. the mapping t 1 — > u°(t,x) is ^i-Holder continuous for < 71 < 

min {E i c Li(pM)i 1 }- 

(ii) For each t £ [0, T) a.s. the mapping x 1 — > u (t,x) is ^-Holder continuous for < 72 < p. 
Proof. 

Proof of (i) Using the semigroup property of the Green kernel G a> s, 

G a ,s(t + h,x - y) = I G a! s(t,x-y- z)G a! s(h,z)dz, 
J«. d 

we have 

u°(t + h, x) - u°(t, x) = [G a>s (t + h,x -y) - G a j{t, x - y)\ u (y)dy 

JR d 

.z) ( / G aS (t,x-y) [u (y-z) -u (y)]dy ) dz. 



/ G aj s(h,. 
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By concavity of the of the mapping x i — > xp/po (since p < po), Lemma |A, II the assumption (H.2) 
and Fubini's theorem we obtain for p £ [2,po] 



E\u°{t + h,x) -u°(t,x)\ P < ( sup E\u (y- z)-u (y)\ Po G atS (h,z)dz 

\jR d y m d , 

< c^G a 4h,z)\zr o d z y° . 



p 
I'd 



Now, due to the scaling property of G c 

G aS (h,z)\z\ ppo dz = I h~^ 1 *G aS (l,h~'* : *^*z)\z\ ppa dz 



= h Ztl^ I GaS ( h y)\y\PPO d y. 

JR d 

Using property (iv) of Lemma ll. II we obtain that 



/ G a 4l,y)\yr°dy < cV / G a ,s(l, yi )\ yi \^dy 



r+oo x ppo 
< c n — dy. 



The last integral is convergent for pp < a - Therefore we have proved the assertion (i). 
Proof of (ii) A change of variable gives immediately 

u°(t,x + z) - u°(t,x) = / [G a ,s{t, x + z - y) - G a ,s{t, x - y)] u (y)dy 

JR d 

= / G a> s(t,x -y) \u {y + z) -u (y)]dy. 

JR d 

Applying again concavity of x i — > x p / po , Lemma I A. II for the probability measure G a g(t, x — y)dy, 
the assumption (H.2) and Fubini's theorem, we obtain for all p E [2,po] 

E\u°{t,x + z) - u°{t,x)\ P < 

p 

< cl [ sup E(\uo(y + z)-u (y)\P°)G a> s(t,x-y)dy ) ' 

\jR d y m d J 

< c^G a 4t, X -y)\zr°d y y° =c\zr. 

We complete the proof of the Lemma by Kolmogorov's criterion. □ 
The next proposition studies the Holder regularity of stochastic integrals with respect to the 
martingale measure M . For a given predictable random field V we set 



U(t,x) 



I I G a>s (t-s,x-y)V(s,y)M(ds,dy). 

JO JR d 
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Proposition 3.2 Assume that sup Q<t<T sup x€R d E(\V(t,x)\ p ) is finite for some p large enough. 
Then under (H") we have 

(i) For each x € M. d a.s. the mapping 1 i — ► U(t, x) is j\~Holder continuous for < 71 < (1 — n)/2. 

(ii) For each t £ [0,T] a.s. the mapping x 1 — ► U(t,x) is f^— Holder continuous for < 72 < 

min{a (l-r?)/2,l/2}. 

Proof. 

Proof of (i). We have 

rt+h r 

U(t + h,x) -U(t,x) = / G a} s(t + h - s, x - y)V(s,y)M(ds,dy) 

Jo Jm d 

- I I G a , s (t - s,x - y)V(s,y)M(ds,dy) 
Jo Jm d 

= [G atS (t + h- s,x-y)-G ai s(t- s,x-y)]V(s,y)M(ds,dy) 

JO JR d 

rt+h r 

/ / G a , s (t + h- s,x-y)V(s,y)M(ds,dy). 

Jt jR d 



+ 

For p > 2, the Burkholder inequality yields 

E\U{t + h,x) - U(t,x) 



p 

p 



< c sup E\V(s,y)\ p {h{h,x)y +(I 2 (M)) f 

(s,y)e[0,T]xR d L 



where 

h(h,x) 

and 



f T ds [ \F [G a , 5 (h + s, x - ■) - G a}S (s, x - •)] (£)| 2 KdO 

Jo JR d 



h{Kx)= f ds f \FG a 4s,x--)(0\ 2 K<%)- 

Jo JR d 

Estimation of Ii(h,x). Set 



^ a ,e(t) := exp{-t^|Ci| ai exp(-^sgn(&))}. (3.15) 



i=l 

By the definition of the Green function G a 5 we can write 

h(h,x) = dr \i(j a ^(r)^ a ^(h) - t/j a ^(r)\ 2 n(d^) 

Jo JR d 
T r 

dr / |^(r)| 2 |^W-l|VK)- 
JR d 

From the proof of Lemma ll,2l we know that \ip a ,^{r)\ 2 < exp(— 2rK,S a (£,)). The mean value theorem 
to the function tp a ^(-) between and h implies that 

\^{h) - 1| < hs Q (o. 
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Moreover, for any (3 g]0, 1[, we have \ip a ^{h) - 1| < 2 1_/3 |^ ajf (/i) - 1|^. 
Therefore 

h(h,x)<cph 2 P [ dr f (S a (0) W exp(-2 K rS a (0)KdO- 



Hence, under (H°) Lemma lA. 31 implies that the right hand side of the above inequality is finite for 
any in ]0, (1 - r,)/2[. 

Estimation of l2(h,x). Fubini's theorem and a change of variable lead to 



I 2 (h,x) = [ \ [ h ds\FG ajS (s,x--)(0\' 

JR d UO 



l^a,^(s)| ds 



< 



exp (— 2nrS a (£,)) dr 



By Lemma lA. 41 the last term is bounded by c(h + h@) for any j3 g]0, 1 — r/[. 
Proof of (ii). Let x G R rf and z belongs to a compact subset K of R rf 



f7(i,x + z) - U(t,x) 



/ / [G aj s(t- s,x + z-y) - Ga tS (t - s,x -y)\ 

JO JR d 

xV(s,y)M(ds,dy). 



For p > 2, Burkholder's inequality yields 

E\U(t,x + z) -U(t,x)\ p < c 



sup E|F( S ,y)| p (J(x,z))2 , 

(s,y)e[0,T]xR<* 



where 



J(x, z) 



I ds I \T[G a j{s,x + z 
Jo Jw. d 



G a>s (s,x--)} (OlV(^) 



The property T (/) (£ + a) = T (e 2l7r { a >')/) (£) together with Lemma TA.5I imply 



J(x, z) 



f ds f |^(s)[exp(2Mr<£,*))-l]|V#) 
JR d 



< 4/ / drexp(-2Kr5 a (0)[l-cos(27r(^,z))]/x(t^) 



R d Jo 



for all < (3 < min{a (l - tj)/2, 1/2}. 



□ 



Proof of Theorem 13.11 Let u be a solution to equation Q1.3JI . For any x G M d , z £ K and 

i G [0,T], and for p G [2,p ], 
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E\u(t + h,x + z)- u(t, x)\ p < cE\u°(t + h,x + z) - u°(t, x)\ p 

+ce( [ [ G a ,s(t + h- s,x + z-y)a(u(s,y))M(ds,dy) 
JO JR d 

- / G a ,s{t - s,x - y)a (u(s,y)) M(ds,dy) 
JO JR d 

rt+h n 

\ \ G a ,s(t + h-s,x + z-y)b(u(s,y))dsdy 

Jo jR d 

- / G at s(t - s,x -y)b(u(s,y))dsdy 
JO JR d 



+cE 



= A 1 + A 2 + A 3 . 

The terms A\, A2 are estimated by Proposition 13.11 and 13.21 More precisely we get 

A 1 + A 2 <c( /i Etl % + h^ p + + \\z\\i 2P 



for /?i < iy2 and p 2 < min ( pi MizZZl^ I). 

Let us now give the estimation of A3. After a change of variable A3 = \B\ P , where 

B = G aj $(h + t - s,x + z - y)b(u(s,y)) dsdy 

Jo JR d 

+ / / G a j(t-s,x-y)[b(u(s + h,y + z))-b(u(s,y))]dsdy. 
Jo JR d 

By Lemma lA.ll to the measure G at $(s, y)dsdy, Lipschitz and linear growth condition of b, we deduce 

A 3 < [ [ G a j(t + h-s,x + z-y)E\b(u(s,y))\ p dsdy 
Jo Jm. d 

+ / G a ,s(t - s,x + z - y)E\b(u(s + h, y + z)) - b (u(s , y))\ p dsdy 
Jo Jsi d 

< c\ ■ h + C2 / sup E \u(s + h, y + z) — u{s, y)\ p ds, 
Jo veR d 



where we have used assertion (i) of Lemma ll. II 

Put ip(s, h, z) = sup vgR d E \u(s + h,y + z) — u(s, y)\ p . Hence, 



sr^d pp ~ ft 

<p(t, h, z) < c 3 (h + h^ i=1 -* + h^ p + + Wzf^) + c 2 / <p(s, h, z)ds. 

Jo 



Therefore by Gronwall's lemma 



sup sup E \u(t + h, x + z) — u(t, x)\ p < 

0<t<T x( zRd 

>d pp 



X^d pp 

< c(h + /i^ =1 * + hP ip + \\z\\Z + 



z 



The conclusion of Theorem 13. II is a consequence of the Kolmogorov continuity criterium. □ 
Remark 3.1 Note that our results of regularity of the solution generalizes those obtained in [26]. 
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4 Smoothness of the law 



We prove that, under non-degeneracy condition on the diffusion coefficient a, the law of u(t, x) 
(solution of Eqg(d, b,a)), has a smooth density for fixed t > 0, x G M. d . 

Since the noise F has a space correlation, the setting of the corresponding stochastic calculus of 
variations is that used in |T7j (see also [T5]). 

Let £ be the inner-product space consisting of functions tp : R d — ► R such that 

r(dx)(v9 * <p)(x) < +00, 

where <p(x) = (p(—x). This space is endowed with the inner product 

( ( P,tp)s= T(dx)(<p * 4>)(x) = / T(dx) dy<p(x-y)ip(-y). 

JR d JR d JR d 

Let TL be the completion of £ with respect the norm derived from (•, Set = L> 2 ([0, T]; Tt), 
this space is a real separable Hilbert space isomorphic to the reproducing kernel Hilbert space of 
the centred Gaussian noise F which can be identified with a Gaussian process {W(h) : h G Ht} as 
follows. For any complete orthonormal system {e^ : k £ N} C £ of the Hilbert space Ti., define 

W k (t)= f I e k (x)F{ds,dx), ken, te[0,T), 
JO JR d 

where the integral must be understood in Dalang's sense. The process {Wfc(i) : t £ [0,T], k G N} 
is a sequence of independent standard Brownian motions, such that for any predictable process X 

[ [ X(s,x)F(ds,dx) = V / (X(s,-),e k (-)) n dW k (s). 

JO JR d k=Q JO 

In particular, for any tp G 5(M rf+1 ) 

^) = E / T (^,-),e fe (-)) w ^ fc ( S ). 
k=o Jo 

For /1 G Hr, set 

W{h) = Y j f (h(s),e k ) H dW k (s). 
k=i Jo 

Thus, we can use the differential Malliavin calculus based on {W(h) : h G TCt} (for more details 
see |19j). Recall that the Sobolev spaces D fc,p are defined by means of iterations of the derivative 
operator D. For a random variable X, D k X defines a Wj! fc -valued random variable if it exists. For 
h G TCt, set D^X = (DX, h)^ and for r G [0, T], D r X defines an element of 7i, which is denoted 
by D r .X. Then for any h G TIt 

D h X= [ (D rt .X,h(r)) n dr. 
Jo 

We write D r)ip X = (D r> .X, p) n for <p G H- 
The main result in this section is 
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Theorem 4.1 Fix t > and x G R d . Assume that the coefficients a and b are C°° with bounded 
derivatives of any order and in addition there exists a > such that o~(z) > a for any z G R. Then 
under (H°) with r\ G (0, |) i/ie Zaw of u(t,x), solution of /ms a density which is infinitely 

differentiable. 

Remark 4.1 If d = 1 and F is a space-time white noise, we obtain existence and smoothness of 
the density of the law for solution to equations studied in J^j and \Tj. In the case where on = 2, 
Si = 0, i = 1, . . . , d, we obtain the result obtained in 

Classical results on the existence and smoothness of the density, using the approach of Malliavin 
calculus, are based on the next proposition. 

Proposition 4.1 W$S[ Let F = (F , . . . ,F m ), m > 1, be a random vector satisfying the following 
conditions: 

(i) F i belongs to D°° = n p >i n fe >i B k ' p for all i = 1, . . . , m. 

(ii) The Malliavin matrix 7^ = ((F l , i ?J )-H)i<i,j<m satisfies 

(det^y 1 G f]LP(Q). 
p>l 

Then, F has an infinitely differentiable density with respect to the Lebesgue measure on R m . 

According to Theorem 2.1 of |15j the proof of Theorem 14. 1\ is a consequence of the following 
Lemma. 

Lemma 4.1 Let t G]0,T]. 

(i) There exists a positive constant c\ such that for any p G [0,t A 1], 

! P ds [ / i(de)iJ r G a ,i( a ,-)(oi 2 >cip < ' 1 

Jo Jm d 

for any 0\ > 1. 

(ii) Assume that (H°) holds for 77 G]0,1[. Then there exists a positive constant C2 such that for 

any p G [0, t A 1], 

j P ds [ v(dO\FG a>6 (s,-)(0\ 2 <c 2P d \ 
Jo JR d 

for any 62 G]0, 1 — rj\. 
Proof. 

(i) Using the lower bound of (|1.4|) . we have, 

pp(d£) 



f ds [ p(dO\FG a 4s,-)(0\ 2 > [ T 

Jo JR d JR d 1 



J R d 1 + b a (t,) 

for any 9\ > 1. 

(ii) Using the upper bound of (|1.4|) . we obtain for any p G [0, t A 1], 

p(d® 



f 'ds [ p(dt)\FG a 4 S ,-)(0\ 2 < 2p f - 
Jo JR d JR d 1 



< cp 1 - 11 



p(dO 



(i+Etil^h)" 

and we complete the proof of (ii) by the assumption (H"). □ 
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A Appendix 

In this last section, we collect some of the technical results which have been used the in previous 
sections. 

Lemma A.l Let f , h be two functions defined on W 1 and \i a positive measure such that f ■ h £ 
L 1 (^i). Then, for all q > 1, we have: 



f.\h\dfi 



< / \ f\ q -\h\dfi) / \h\dfi 



9-1 



Proof. Set v = \h\dfi, then the result follows from the Holder inequality applied to J fdv. □ 
We state below without proof a version of GronwalFs Lemma that plays a crucial role in the 
proof of the existence and uniqueness results. 

Lemma A. 2 (Jj^j, Lemma 15) Let g : [0, T] — ► R + be a non-negative function such that g(s)ds < 
+oo. Then, there is a sequence {a n ,n £ N} of non-negative real numbers such that for all p > 1, 



2~2n=i ( a n) 1 ^ P < and having the following property: Let (f n ,n £ N) be a sequence of 

negative functions on [0, T], k\, k 2 be non-negative numbers such that for < t < T, 

fn{t)<h+ [ (k 2 + fn-l(s))g{t-s)ds. 

Jo 

If snp 0<s<T f (s) = c, then forn>\, 



non- 



n-l 



fn{t) <h + (ki + k 2 ) ^2 Oi + (k 2 + c)a n . 



i=l 



In particular, sup n>0 sup 0<t<T f n {t) < +oo. 

If k\ = k 2 = 0, then X] n >o (/n(*)) converges uniformly on [0,T]. 
Lemma A. 3 Let n and j3 be in (0, 1), if (H^) holds then 

[ T dr [ exp (-2r K S a (0) (S Q (0) 2/ V(<£) 

JO JR d 

is finite for all (3 in ]0, (1 — rj)/2[. 

Remark A.l Since the spectral measure fi is non-trivial positive tempered measure, we can ensure 
that there exist positive constants c\, c 2 and k such that 

ci < fj,(d£) < c 2 

J\G\<k 

Proof. Let (3 e]0,l[, set 

I(h,x) = Fdrf exp(-2 K rS a (0)(S a (Z)) 2/3 

Jo JR d 
= h + h 
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where 



and 



'1 



dr 



e W (-2KrS a (0)(S a (0) 2 ^(dO 



{s Q (0>i} 



ds exp (— 2KrS a (^)) 



f 

J {S a (€)>!} 



1 - exp(-2r«g a (g)) 

2~K 



'{s Q (5)>i} 

Clearly 7i is finite. And the term I 2 is bounded by 

1 



K J {S a (.H)>l} i 1 + S a{0) 



<- sup (S a (0) 2 '- 1+V f TTT^W- 

Now, choose /? in ]0, (1 — 77) /2[ and by the hypothesis (H^) the last term is finite. 
Lemma A. 4 If (H^) ZioWs i/ien 

rt+h r 

sup / / \TG a ^{t + h- s,x - -){i)\ 2 ii{di) <c{h + h p ) 

0<t<T Jt JR d 

for any (5 G]0, 1 — rj[. 

Proof. Let h be in [0, T] such that t + h belongs to be in [0, T] and set 

I 2 (h)= f ds [ \FG atS (s,x--)(0\ 2 »{dO. 

JO JR d 

By the definition of G a ,s we recall that JG Q] j(s,s — •)(£) = ip a ^(s) where ip a ^(s) is given by 
(|3.15|) . Using Fubini's theorem we can write 



h{h) 



\lp a As)\ ds 



M(de) 



< 



exp (— 2rer5' a (£)) dr 
h.i(h)+I 2 . 2 (h), 



where 



and 



h.i(h) := [ \ [ h exp(-2KrS a (0)dr 

J{s a (i)<i} i J o 



/ [/ 



{s a (t)>i} 



exp (— 2nrS a (£ i )) dr 

1 — exp(— 2nhS a (^)) 
2kS q (0 
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It is easy to see that h.i(h) is bounded by c • h. 
Moreover for each (3 e]0, 1[ we have 

|1 - exp (-2KhS a (g))\ < 2 l -P |1 - exp (-2KhS a (£))f . 

Using the mean value theorem we get 

J{s a (0>i} 

Choosing (3 G]0, 1 — n[ and use hypothesis (H") we show easily that 

/ (s Q (0) /3_1 M^)<+oo. 

J{s a (0>i} 

Lemma A. 5 Let z belongs to a compact subset K ofW 1 . Under (H") we have 
[ [ dr exp (-2 K rS a (0)[^-- cos (27r{^z))]^)<c(\\z\\ 00 + \\z 

JR d JO 

for any < (3 < min{l/2, a (l - rj)/2}. 
Proof. Let us put 

J(z) = J^(^£ drexp(-2KrS a (0)d S ^[l-cos(27T^,z))}fi(dO 
= Ji(z) + J 2 (z), 



2/3 \ 
oo / 



where 



Ji(z) : =/ dr f exp(-2Kr5 Q (0)[l-cos(27r^,z))]M(de) 
J J{S a (£)<l} 

< [ [1- cos (2tt<£, 

J\S„ (£)<!} 



and 



■M* : =/ P 9 \ q m { " 1-cos 2tt g,z M dfl 



'{Sa(0>l} 

< 2 1 - 2 / 3 /■ |l-cos(27r(e,z))| 



for any [3 €]0, l/2[. 

By the mean value theorem 



Ji(z) < c |zd < c ■ d ■ \\z\ 



i=l 
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and 

Hz) < c P , K ! (E?=1 'y )2 % (de) 

<r IMI 2 / 3 f (Ei=i I6I) 2/3 

^{S«(0>1} 

If 1 < «j < 2 for i = 1, . . . , io, set ao,i = m i n i<i<io{ a «} > !• Hence 

2/3 2/3 

(*o \ M / io \ «o,i / l o \ «o,i 

Ei&i) <^(Ei^r M ) <% d (Ei&H > (A.i6) 

where we have used the convexity property of function x i — ► x 7 for 7 > 1 and the inequality 

If < aj < 1 for i = zq, . . . , d, set «o,2 = rnnij <j<d{aj}. Use the fact that Ix) 1 /" 4 < |x| 1//qo ' 2 for 
|x| > 1 and the super-additivity of the function x 1 — ► x 7 for 7 > 1 to obtain 

/ d \ 2 P / d _^\ 2/3 / d \^ 

\i=i / \i=io / \i=io / 

Recall that «o = ^^i<i<d{ a i}y then qo = min{ao,i, 00,2}- Combining l)A.16|) and ()A.17|) with 
concavity of the mapping x 1 — ► x 7 (7 < 1) for x > 1, we end up with 

/ d \ 2/3 / d \^ 



for 2(3 < ceo. Consequently 



/ 



E&i £i6r 

vi=0 / Vi=0 



(E-=i ie 



^{SaK)>l} 

which is finite according to (H") for all < j3 < ao(l — rj)/2. □ 
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